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ABSTRACT 

Hot accretion flows such as advection-dominated accretion flows are generally optically thin 
' in the radial direction. Thus photons generated at some radii can cool or heat electrons at other 

radii via Compton scattering. Such global Compton scattering has previously been shown to 
be important for the dynamics of accretion flows. Here, we extend previous treatments of this 
problem by using accurate global general relativistic Monte Carlo simulations. We focus on 
an inner region of the accretion flow (R < 600/? g ), for which we obtain a global self-consistent 
solution. As compared to the initial, not self-consistent solution, the final solution has both 
the cooling rate and the electron temperature significantly reduced at radii > 10 gravitational 
radii. On the other hand, the radiation spectrum of the self-consistent solution has the shape 
similar to that of the initial iteration, except for the high-energy cut-off being at an energy 
lower by a factor of ~ 2 and the bolometric luminosity decreased by a factor of ~ 2. We 
also compare the global Compton scattering model with local models in spherical and slab 
geometry. We find that the slab model approximates the global model significantly better than 
the spherical one. Still, neither local model gives a good approximation to the radial profile 
of the cooling rate, and the differences can be up to two orders of magnitude. The local slab 
model underestimates the cooling rate at outer regions whereas it overestimates that rate at 
inner regions. We compare our modelling results to observed hard-state spectra of black-hole 
binaries and find an overall good agreement provided any disc outflow is weak. We find that 
general-relativistic effects in flows which dynamics is modified by global Comptonization is 
crucial in approaching this agreement. 
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1 INTRODUCTION 

The lowest-order interaction between photons and electrons, 
namely Compton scattering, results in momentum and energy ex- 
change. It is a common and important process in astrophysics. The 
momentum exchange regulates accreting gases to accrete only be- 
low a critical value, corresponding to the Eddington limit on the 
luminosity, L E (for spherical accretion). The importance of energy 
exchange manifests itself in two aspects. For photons, Compton 
up-scattering by energetic electrons is the main mechanism to pro- 
duce X-ray emission in various astrophysical systems, and Comp- 
ton down-scattering by low-energy electrons often leads to mod- 
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ification of X-ray spectra. For electrons, they either gain or lose 
energy, depending on the average photon energy. This process is 
important, in particular, in hot accretion flows. 

Hot accretion flows such as advection-dominated accretion 
flows (ADAF, e.g., Narayan & Yi 1994; Abramowicz et al. 1995) 
are optically thin in both vertical and radial directions. Thus, a pho- 
ton can travel a long distance before being absorbed or scattered 
(e.g., Narayan & Yi 1994; Narayan, Mahadevan & Quataert 1998). 
Then, photons produced in one location can either cool or heat the 
flow in another location, and Compton effect is not local but global. 
The extension of ADAF to higher accretion rates (but still below 
the Eddington limit), namely luminous hot accretion flows (LHAF, 
Yuan 2001 hereafter Y01; Yuan 2003), should be even more af- 
fected by the global Compton scattering, since its optical depths are 
higher and the radiation fields are stronger than those of an ADAF. 
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Global Compton cooling effects due to external photon 
sources like the surface of a neutron star have been investigated 
in the past (e.g., Narayan & Yi 1995). In this study, we concentrate 
on pure hot accretion flows (ADAF or LHAF) where all the seed 
photons for Compton scattering processes are generated within the 
hot plasma. We note that if the geometry of an accretion flow is an 
inner hot flow plus an outer truncated thin disc (e.g., Esin, McClin- 
tock & Narayan 1997), cooling by the outer thin disc (Shakura & 
Sunyaev 1973) should also be included; however, we neglect it here 
for simplicity. 

While Comptonization in the vertical direction is included 
in the standard ADAF/LHAF models by a local one-zone treat- 
ment (e.g., Narayan & Yi 1995; Manmoto, Mineshige & Kusunose 
1997), any non-local radiative transfer in the radial direction, es- 
pecially Comptonization, has been neglected in almost all of the 
previous works in this field. Global Comptonization in hot accre- 
tion flows has been considered, as far as we know, only by Esin 
(1997), Kurpiewski & Jaroszyhski (1999), Park & Ostriker (1999, 
2001, 2007), and recently by Yuan, Xie & Ostriker (2009, hereafter 
YXO09). We note that only YXO09 and Kurpiewski & Jaroszyfiski 
(1999) attempted to find a global solution to this problem rather 
than to use a self-similar approach, which fails at the inner regions, 
where most of the radiation is generated. YXO09 used an iteration 
method to find, for the first time, the radiative cooling rate mu- 
tually consistent between the dynamics and radiation of the flow. 
They have found that Compton heating and cooling dominates out- 
side and inside, respectively, of ~ 10 4 /? g , where R g = GM BH /c 2 
is the gravitational radius and Mbh is the black-hole mass. They 
have also found that the relative effect of global Comptonization 
increases with the increasing accretion rate. 

We note, however, that YXO09 have concentrated on the 
Compton heating of the outer flow. Also, the only photon produc- 
tion rate they used was from the intrinsic emission (synchrotron 
and bremsstrahlung) and their local Comptonization at a given ra- 
dius. Thus, although they did calculate the effect of the Compton 
energy exchange (using the exact relativistic treatment of Guilbert 
1986) on the electron temperature in a given radius including pho- 
tons produced at all radii, they did not take into account the photons 
generated after the non-local scattering. Thus, non-locally scattered 
photons were assumed to leave the flow. Since outer parts of hot 
flows are very optically thin, this is a reasonable assumption for 
calculating their Compton heating. 

On the other hand, here we consider Compton cooling of an 
inner part of a hot flow, where the Thomson optical depth becomes 
comparable to unity at moderate accretion rates. Thus, we cannot 
use the method of YXO09. Instead, we use a Monte Carlo (here- 
after MC) method for Comptonization in the hot flow. Our method 
is directly based on that of Niedzwiecki (2005) and Niedzwiecki & 
Zdziarski (2006, hereafter NZ06), which is a generalization of the 
method of Pozdnyakov, Sobol' & Sunyaev (1983) and Gorecki & 
Wilczewski (1984) to include effects of special and general relativ- 
ity (hereafter GR) and the bulk motion of the flow. Comptoniza- 
tion by bulk motion (considering inflow only) was earlier stud- 
ied by, e.g., Titarchuk, Mastichiadis & Kylafis (1997) and Laurent 
& Titarchuk (1999; see also NZ06 for a critical discussion). We 
assume that electrons have a relativistic Maxwellian distribution, 
without any non-thermal contribution. In our dynamical treatment, 
we use a one-dimensional model, which assumes the flow to be 
both azimuthally and tangentially uniform, with the scale height 
given by hydrostatic equilibrium. 

MC simulations for Compton scattering process in a hot ac- 
cretion flow in the Kerr metric was first done by Kurpiewski & 



Jaroszyfiski (1999). However, their dynamical model was based on 
the fully advection-dominated flow (i.e., with negligible radiative 
cooling), which is clearly not the case for flows with higher accre- 
tion rates, which we investigate here. 

We describe our treatment of the problem in Section [2] We 
present a comparison of various Comptonization models in Section 
[3] In Section[4] we discuss the effects of outflows and viscous elec- 
tron heating on the accretion solutions. The self-consistent solution 
for our chosen accretion rate is discussed in Section[5] Our results 
are compared with observations in Section [6] and conclusions are 
given in Section|7] 



2 THE MODEL 

Here, we state the problem we are solving, and present the dy- 
namical equations and the Compton scattering method we use. 
We focus on a hot accretion flow around a stellar black hole with 
Mbh = IOMq. The outer boundary is set at R out = 600R g and the 
accretion rate there is Mq. 

In order to study the effect of global Comptonization, we need 
to know the structure of the flow, i.e., the density, temperature and 
velocity as functions of the radius. We iterate between the dynam- 
ical solutions and the MC Comptonization results to obtain the 
self-consistent Compton cooling rate. We start from solving the dy- 
namical structure without considering the global Comptonization, 
which is customary in most of the previous work on hot accretion 
flows (e.g., Narayan & Yi 1995; Manmoto et al. 1997; Quataert & 
Narayan 1999; Yuan, Quataert & Narayan 2003). Then we calculate 
the radially-dependent Compton cooling by the MC simulations. 
We then use this rate to calculate again the dynamical structure. We 
repeat the above procedure until it converges. 



2.1 Dynamical equations of hot accretion flows 

The scale height, H, of the hot accretion flow in our solution de- 
pends on the radius. It is given by hydrostatic equilibrium, 
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where p is the total pressure, which includes contributions from 
electrons, ions and magnetic field (see below), and p is the mass 
density, where we assume cosmic abundances with the H mass frac- 
tion of X = 0.75. We assume that the density is tangentially uniform 
from the disc mid-plane up to ±H . We adopt the pseudo-Newtonian 
potential of Paczyhski & Wiita (1980), in which Sl K (the Keplerian 
angular velocity) is given by, 
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Because of the effect of outflow/convection, the mass inflow 
rate is function of radius and we assume it to be given by 
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where v (< for the inflow) is the radial velocity. Hereafter, we 
adopt the unit of Me = Le/c 2 for Mq, where Le = 47tGMbh'"pC/0't 
is the Eddington luminosity for pure H, m p is the proton mass and 
cr x is the Thomson cross section. The value of s is well constrained 
in the case of the supermassive black hole in our Galactic Cen- 
tre (Yuan, Quataert & Narayan 2003), and we set s = 0.3 follow- 
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ing that work. We discuss the validity and consequences of this 
assumption in Sections |4]and[6]below. 
The energy equation for electrons is, 
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where £ c is the specific internal energy of electrons, p c is the elec- 
tron pressure, and gi e is the electron heating rate per unit volume 
by ions via Coulomb collisions. The electrons can also be directly 
heated at the rate 6q v i s , where q vls is the viscous heating rate per unit 
volume. We assume 6 = 0.5, again following the detailed study of 
the Galactic Centre accretion flow of Yuan et al. (2003). The elec- 
trons are cooled by synchrotron and bremsstrahlung, q syn+ bi, and by 
the Compton scattering process, q comp . The latter is the net Comp- 
ton cooling rate, which is, by definition, 



E) dE. 



(5) 



Here, J E is the mean intensity per unit photon energy, E, in the 
local rest frame at R, cr(E, T s ) is the Compton cross section aver- 
aged over the electron Maxwellian distribution at the temperature, 
T e , and (E t ) is the average photon energy (in the same frame) af- 
ter scattering. However, since we use MC simulations (see Section 
12.21 ) to find our self-consistent solution, we calculate the (equiva- 
lent) Compton energy exchange rate directly from the simulations 
instead of using equation (O. This rate differs from the local Comp- 
ton cooling rate (commonly used in other studies of hot accretion 
flows) in that J E includes now photons either produced or scattered 
at all radii that reach a given point at R. It can also become negative, 
i.e., corresponding to the radiative heating by photons, especially in 
outer regions of the flow ( > 10 4 i? g ; e.g., YXO09). However, since 
we concentrate here on inner flows, we find q comp > in all cases 
considered here. 

We assume the ratio of the gas pressure (electron and ion) 
to the magnetic pressure, /?b = 9, i.e., the magnetic pressure of a 
l/10th of the total pressure. This parameter determines the strength 
of the magnetic field in the accretion flow, and the synchrotron 
emission can then be determined. For it, we follow the method of 
Narayan & Yi (1995) and Manmoto et al. (1997). 

For completeness, we list below the other dynamical equa- 
tions, i.e., equations for the radial and angular momenta and for the 
ion energy (see, e.g., Yuan et al. 2003; Manmoto et al. 1997), 

dv , , 1 dp 

v — = -SliR + a~R -, (6) 

dR K pdR 

v(OR 2 - f) = -aRp/p, (7) 

where Q is the angular velocity of the flow, j is the specific angular 
momentum of the accretion flow when it crosses the horizon, and 
a is the viscous parameter, which we set here as a = 0.3. 

We note that except for using the pseudo-Newtonian poten- 
tial, our dynamical equations are non-relativistic (NR). The radial 
velocity then un-physically exceeds the speed of light in the vicin- 
ity of the black hole. We have tried two kinds of corrections for 
this. First, in modelling the Galactic Centre, Yuan, Shen & Huang 
(2006) divide the velocity by factor 0.93 exp(4.26R g /R), which 
comes from fitting to corresponding relativistic results. Second, we 
note that a first-order dynamical relativistic correction is to replace 
the NR momentum, mv, by ymv, where y is the Lorentz factor of the 
flow. We then approximate the actual velocity using the assumption 
that the velocities in the NR calculations are the actual ones multi- 



plied by factor y. We have compared these two methods and found 
no significant differences. Here, we use the latter treatment. 

2.2 Monte Carlo method for Compton scattering 

The MC method used here is fully GR in the Kerr metric, and 
it is described in Niedzwiecki (2005) and NZ06. However, in the 
present paper we consider only a non-rotating black hole. We thus 
use the following dimensionless parameters, 



R 

— , I'- 
ll. 



Lc 

E m fR„ 



EjrfRg 



c 



p* ■ 



(9) 



where L = R x p is the vector of the photon angular momentum, 
P = (Pr, Pe, Pip) is the photon momentum vector, R = (R, 0, 0), is 
the component of L parallel to the symmetry axis, i.e., L$ = L cos 
and f is the (dihedral) angle between the plane of photon motion 
and the equatorial plane. Then, E- mf is the photon energy at infinity, 
v' and v* are the radial and azimuthal velocities of the flow, respec- 
tively, with respect to local static observers (i.e., observers with r, 
6 and = constant). We use the relativity-corrected (as described 
above) values of v and RQ. as v' and v*, respectively. 

The tangential motion of the accretion flow is neglected, v 9 = 
0. The tangential (or vertical) position of a photon within the flow 
is given by 6. The energy generation and exchange rates, q, are 
defined in the rest frame of the flow. Seed photons are randomly 
generated following the radial distribution of q syu+ bi(r) and are as- 
sumed to be isotropic in the local rest frame. The flow is assumed to 
be uniform tangentially, and thus photons are generated uniformly 
in 9. Then, we follow the photon trajectory in the curved space- 
time through subsequent scatterings until the photon either crosses 
the event horizon or escapes the flow. In the latter case, it then con- 
tributes to the observed angle-dependent spectrum. 

Note that the optical depth in a relativistic flow around a black 
hole is affected by both special and general relativistic effects (see 
NZ06 for a detailed discussion). The following form is valid for 
electron scattering in the Schwarzschild metric, in a spiral rotating 
flow with v e = 0, 
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where f is the affine parameter. The upper and lower sign above 
corresponds to ingoing (i.e. with dr/d£ < 0) and outgoing photons, 
respectively (note that p < 0). 

The simulation of Compton scattering is performed in the lo- 
cal rest frame using the MC method of Gorecki & Wilczewski 
(1984; see also Pozdnyakov et al. 1983). The electron temperature, 
the density, and the velocity vector are found in each iteration from 
the hydrodynamical solution by interpolating between values cal- 
culated at 57 logarithmically spaced radii, r, (such a grid is suffi- 
cient given a smooth distribution of the flow parameters). 

In order to allow us to compare our results with some previ- 
ous models, which usually neglect general and special relativistic 
effects and use local approximations for Comptonization, we per- 
form also calculations for local models. We define a set of 57 uni- 
form, independent spheres or infinite slabs with the n c , T s , and the 
seed spectrum (synchrotron and bremsstrahlung) given for each r,. 
The sphere radius or the slab half-thickness are given by the corre- 
sponding scale height of the flow, Then, the Thomson optical 
depth corresponding to the sphere radius or the slab half depth is 
tjj = <T7n c (rj)Hj. We then use the MC Comptonization method 
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Figure 1. The cooling rates times volume, cfR 2 H, for models based on the 
dynamical properties of the initial solution with Mo = 0.5Me. The black 
dashed curve gives the rate due to the synchrotron and bremsstrahlung emis- 
sion. The thick magenta and black solid curves correspond to the local slab 
model and global GR Compton scattering model, respectively, both of them 
obtained using the MC method. The thin red solid curve gives the cool- 
ing rate using the local slab model of E96 but with the optical depth from 
D91 for their 1-eV seed photon case. The blue short-long dashes and green 
dashes correspond to the local approximation with the model of D9 1 for a 
sphere and slab, respectively, with the coefficients tabulated for 1-eV seed 
photons. See Section[5]for details. 



in flat space for each slab or sphere. The local Compton cooling 
rate can then be determined for each slab or sphere, while the to- 
tal luminosity is the sum of the radiation from all the 57 slabs or 
spheres. 



3 COMPARISON OF COMPTONIZATION MODELS OF 
HOT FLOWS 

We first consider some Comptonization models used in other works 
and compare them with our results. A quantity often used to quan- 
tify the effect of thermal Compton up-scattering is the amplification 
factor, 77, by which the power in initial seed photons is amplified 
through Compton scattering. In our formalism, 



rj(r) ■ 



?syn+br('") ~^ ^comp('") 



(ID 



In local models, it is assumed that only photons produced at r are 
Comptonized at this radius. The amplification factor of the Comp- 
tonization depends on the local electron temperature, T e , the Thom- 
son optical depth, t t , and on the initial photon energy, E\. Here, we 
use the formalisms of Esin et al. (1996, hereafter E96) and of Der- 
mer, Liang & Canfield (1991, hereafter D91), averaging their in- 
dependent 77 over our seed spectra to get the corresponding Comp- 
ton cooling rate. We note whereas E96 give an approximate formula 
that is supposed to be valid for any E,, equations (1), (20), and ta- 
ble 1 of D91, which we use here, contain coefficients tabulated at 
E\ = 1 eV and 1 keV only. Thus, their formula is accurate only for 
those values of E\ in spite of the appearance of E\ as an argument 
of their equation (1). 



Fig.Q]shows the cooling rates for the dynamical solution given 
by the initial iteration in Section|5]at M = 0.5M E . In this initial 
iteration only, local slab Comptonization is applied (in order to get 
the first iteration of the dynamical solution), and the Comptoniza- 
tion cooling (shown as thin red solid curve) is treated using equa- 
tion (A10) of E96 and equation (20) of D91, with coefficients taken 
from the E\ = 1 eV slab case. We use here the 1 eV case since most 
of our synchrotron seed photons are close to this energy. (The re- 
sults using the original formulae of E96 are very similar to the thin 
solid curve and thus not shown here.) In addition to Comptoniza- 
tion cooling rates, we also show the cooling rate due to the intrinsic 
seed emission via synchrotron and bremsstrahlung processes. Then, 
the ratio of Comptonization cooling rate to the seed one equals to 
r] — 1. We also show here the analytical models of D91 for sphere 
and slab, and compare them to our local NR results and the global 
GR MC results. 

The thick black solid curve shows global GR MC results based 
on the initial dynamical iteration. We can see that the worst model 
compared to these accurate results is the spherical one (blue short- 
long dashes). This is because a sphere is a very poor approximation 
to the geometry of an accretion flow, more similar to a slab. A key 
factor affecting Comptonization in hot accretion flows is that the 
optical depth in the radial direction is much larger than that in the 
vertical direction (which is ^ 0.7 from horizon to the outer bound- 
ary for the initial solution). Only models with a slab geometry can 
account for this property. The assumption of a spherical geometry, 
which ignores the large difference in the radial and vertical optical 
depths, leads to significant underestimation of the probability of 
scattering, and thus under-prediction of both the magnitude of the 
Compton cooling rate by about an order of magnitude (see Fig. [I}, 
and the hardness of Comptonization spectrum, see Section[5]below. 

As we see in Fig. [T] although the local slab model provides 
much better approximation to the radial profile of the cooling rate 
than the spherical model, deviations with respect to the global GR 
calculations remain highly significant. The local slab approxima- 
tion, in either version of E96 or D91, has several deficiencies, 
which we discuss below. 

First, the seed photons in the formalism of D91 are assumed 
to be at either 1 eV or 1 keV. This is obviously not the case in gen- 
eral. Especially when most of the power in seed photons is at high 
energies, as it is the case for bremsstrahlung, that formula becomes 
highly inaccurate. Since typical electron temperatures in hot accre- 
tion flows are mildly relativistic, scattering is affected by Klein- 
Nishina effects. Those effects are also not taken into account by 
E96. Thus the formulae for the amplification factors from those 
works usually lead to an overestimation of the Compton cooling 
rate when bremsstrahlung radiation is the main seed photon source. 
This can be seen in Fig. Q] at R > 200R g , where the local analytical 
results are way above the local slab MC results. 

Second, even more importantly, Comptonization is not local. 
As shown in Fig. [T] Comptonization can increase the cooling rate 
by up to two orders of magnitude, i.e., the process is very important. 
In particular, we see that the cooling rate for global Comptoniza- 
tion is much less centrally concentrated than that for the local treat- 
ments. The seed photon production rate itself is centrally concen- 
trated; however, those photons travel to outer regions, increasing 
the cooling rate there. The local slab model neglects this effect. 

Third, strong gravity leads to capture of photons generated in 
the innermost region of the accretion flow, and redshifts escaping 
photons. This will be further discussed in Section[5]below. 

Fourth, Comptonization on the flow bulk motion may, in prin- 
ciple, be important in addition to Comptonization by electron ran- 
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dom motion. Here, we have both rotation and inflow. The rotation 
remains slightly sub-Keplerian throughout most of the flow, and 
the radial velocity, v aSl K R (e.g., Narayan & Yi 1994) is signif- 
icantly smaller than the free-fall velocity. Then, both components 
of velocity become relativistic only in the innermost part, within 
a few R g . The bulk plasma motion can affect the spectrum of the 
observed radiation through the following effects, (i) Collimation 
of the radiation along the direction of plasma motion, which in a 
converging flow around a black hole results in the increase of the 
number of photons captured by the black hole. On the other hand, 
photons emitted outward have an increased scattering probability 
due to the inflow. This increases photon trapping in optically-thick 
flows. Even in optically thin flows, photons after those scatterings 
are collimated inward, and mostly captured, (ii) The same effect 
changes the photon energy after scattering, and the change now 
is due to both the bulk and thermal motions of the electrons. A 
condition for the bulk motion to be an important source of energy 
can be written as (y - l);n c c 2 /(3kT e ) > 1 (e.g., Blandford & Payne 
1981), which is virtually never satisfied in advection-dominated ac- 
cretion flows (see Jaroszynski 2001). In our simulations, we have 
compared the power transferred to photons in the plasma rest frame 
(i.e., from thermal motion only) and in the local static frame (with a 
contribution of bulk motion), and have found that the latter is larger 
by only ~ 10 per cent greater around lOR g . Only at = 23R g the flow 
velocity becomes highly relativistic (and thus motion contributes 
appreciably), but the contribution of this region to the total radiated 
power is tiny, as seen in Fig.[T] At > 10i? g , the effect of bulk motion 
is completely negligible. Overall, we conclude that the bulk motion 
Comptonization has only a minor effect for hot accretion flows. 



4 EFFECT OF OUTFLOWS AND VISCOUS HEATING OF 
ELECTRONS ON ADAF SOLUTIONS 

Taking into account of disc outflows (e.g., Stone, Pringle & Begel- 
man 1999; Blandford & Begelman 1999) and of viscous heating 
of electrons (e.g., Bisnovatyi-Kogan & Lovelace 1997; Quataert & 
Gruzinov 1999; Sharma et al. 2007) represent two progresses of 
the hot accretion theory since its initial formulation. Including of 
either of these possible effects influences accretion solutions in a 
major way. Since in this work we investigate how global Comp- 
tonization affects those solutions, for completeness and to enable 
comparison among all of these effects, we also show how outflows 
and viscous electron heating modify the accretion flow structure. 

Fig.[2]shows a comparison between the ADAF model without 
an outflow, s = 0, and with no viscous electron heating, S = 0, and 
models with either s or 6 or both being > 0. These solutions are ob- 
tained with local Comptonization using the analytic approximation 
of E96 and D91 used by us (see Section[3]l. We see first that adding 
viscous electron heating has a relatively minor effect. It somewhat 
increases t t and T e and slightly decreases T,. On the other hand, 
including outflows leads to a dramatic decrease of tt and relatively 
strong increase and decrease of T e and Tj, respectively, with the ef- 
fects being much stronger in inner flow regions. Then, setting both 
s = 0.3, 6 = 0.5 results in a superposition of the individual effects. 

These effects can be understood as follows. Increasing 6 in- 
creases and decreases the power supply to electrons and ions, re- 
spectively. An increase of s obviously strongly reduces the density 
and tt. This reduces the radiative cooling (roughly oc p 2 ) faster than 
the viscous heating (oc p), which results in an increase of T e . The 
corresponding decrease of Tj is due to the reduced compression 
work caused by the shallower density profile. 



5 RESULTS 

Fig. [3] presents comparison between the initial (short-dashed 
curves) and the final self-consistent (solid curves) global GR MC 
solutions at Mo = 0.5Me. Fig.(3{a) shows the electron and ion, T[, 
temperatures. We see that compared to the initial solution, T c de- 
creases significantly at > 10/?„, up to a factor of ~ 2 at ( 100-200)^ 5 . 
This is an effect of the increased cooling in that region, as we see 
in Fig.[3]c). On the other hand, there is almost no change in T\. 

The height-to-radius ratio is shown in Fig. [3jb). We find that 
the height of the flow is nearly unchanged at <20i? g , and slightly 
decreases at larger R, which is due to a decrease of T- l (and T s ) 
there. We also find only relatively small changes in v and n c . The 
vertical Thomson optical depth shown in the inset of this panel in- 
creases slightly, while it remains nearly unchanged at < WR g . The 
horizontal optical depth of the final solution is 0.75. 

Fig. [3]c) compares the cooling rates for the global GR model 
and also for the local slab model. The final solutions differ signif- 
icantly from the initial ones. Also, the global MC solutions differ 
very strongly from those for the local slab model. The final self- 
consistent profile is relatively flat, namely, drops by only a factor 
of a few from the maximum at ~ 10i? s to 200R g . The flat profile is 
partly due to the global Compton scattering. Also, there is a very 
strong outflow in our model, see equation This strongly reduces 
both M and the cooling rate towards the black hole. 

Fig. |3jd) shows the spectra of the initial and final solutions 
from MC simulations with both the local slab and global GR mod- 
els. We find that the shape of the spectrum remains very similar for 
the initial and final solutions, but the normalization of the spectra of 
the final solutions is about a half of that that of the initial solution. 
This is roughly consistent with the result of YXO09. The bolomet- 
ric luminosity of the final model is L ^ 8.4 x 10 36 erg s -1 , which 
is 6.7 x 10~ 3 Le. If we define the radiative efficiency as L/(MoC 2 ), it 
will be 0.013. However, we stress that most of the potential energy 
of the flow at R aul is lost to the assumed outflow. 

For either the initial or final solution, the spectral differences 
between the global GR and local slab models are relatively small, 
in particular much smaller than the differences in the cooling rate. 
This is because each spectrum is the sum of the radiation pro- 
duced at all radii, which averages the local differences. For the final 
self-consistent solution, the difference between the global and lo- 
cal models is significant mostly at the high energy end, > 100 keV. 
This is because this part of the spectrum is produced by Comp- 
tonization in the innermost part of the flow, where the capture by 
the black hole and other GR effects are strong, which reduces the 
energy of the high-energy cut-off. A corresponding difference for 
the initial solutions does not appear because the high-energy radi- 
ation does not originate from the innermost region, as we can see 
in Fig.[3]c). The differences in the soft X-ray/UV bands are also 
caused by the capture and redshift effects. We also see that the lo- 
cal spherical MC model leads to a strong under-prediction of the 
hardness of the spectrum, see the dotted curve in Fig.[3jd). 

When the accretion rate increases (but still L <k Le), the radia- 
tive cooling will exceed the local viscous heating, and the accretion 
flow will enter into the LHAF regime (Y01). The accretion flow can 
still remain hot because of the heating by the compression work 
(Y01) and it is the sum of the viscous dissipation and the compres- 
sion work that balances the radiative cooling. Above a certain M a \u 
the radiative cooling becomes so strong that the flow cannot remain 
entirely hot. Such a flow below a certain R (see, e.g., the long dashes 
in fig. 1 of Y01) will either collapse into a thin disc or switch into a 
two-phase flow, with cold clumps embedded in the hot phase (Y01, 
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Figure 2. The effect of an outflow (measured by s) and viscous heating of electrons (at the fractional rate of S) on the ADAF solutions. The dot-dashes, short 
dashes, long dashes and solid curves show the original ADAF model (s = S = 0), and models with strong viscous electron heating and no outflow (s = 0, 
6 = 0.5), with no viscous electron heating but strong outflow (s = 0.3, <S = 0), with both phenomena included (s = 0.3, <S = 0.5). In all models, Comptonization 
is treated analytically (D91; E96), and Mq = 0.5Me- (a) The radial profiles of the electron (lower curves) and ion (upper curves) temperature, (b) The radial 
profile of the vertical optical depth. 




Figure 3. Comparison between the initial (short dashes) and final self-consistent (solid curves) MC solutions with Mo = 0.5Mg. (a) The radial profiles of the 
electron (lower curves) and ion (upper curves) temperatures, (b) The ratio of the scale height to the radius, and the vertical optical depth shown in the inset, (c) 
The radiative cooling rate, (d) The spectrum. In (c-d), the (red) long dashed and dot-dashed curves correspond to local-slab MC simulations of the initial and 
final solutions, respectively. The blue dotted curve in (d) shows the spectrum from the local-sphere MC simulations for the dynamical structure of the initial 
solution. [Note that the short and long dashes in (c) are, respectively, the same as the thick solid black and magenta curves in Fig.[T|] 
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Figure 4. Comparison between the radial profiles of the electron (lower 
curves) and ion (upper curves) temperatures for MC calculations for the ini- 
tial (dashed curves) and a solution including global Comptonization (solid 
curves) for Mq = I.OMe- We see there is no self-consistent solution for 
i?< 10 2 i? g due to the very strong cooling. At lower radii, either the disc 
collapses or a two-phase cold/hot flow forms. 



Yuan 2003). The value of M crit decreases when the global Comp- 
ton scattering effect is included. Here, we find that at Mo = I.OMe, 
the flow collapses at < 100W g , as shown in Fig. [4] Thus, the criti- 
cal accretion rate for our assumed flow with the strong outflow at 
s = 0.3 is given by 0.5 < M clit /M E < 1.0. This M clit will decrease 
even further if the cooling due to the photons from outer thin disc 
is included. We also find that M clit is higher in the local-Compton 
model with the cooling rate of E96 and D91 (as used by us above), 
- 1.5Me (which solution is not presented here) than in the global- 
Compton model. This value is then lower than that found by Y01, 
which is in the range of ~ (3-5)Me- (Note that Y01 used the defini- 
tion of Me including an efficiency factor of 0.1.) The discrepancy 
is because both the strong outflow, which suppresses the compres- 
sion work (crucial for LHAF), and the strong direct electron heating 
(6 = 0.5) are included in the present model while they are not in 
Y01. 



6 COMPARISON WITH OBSERVATIONS 

Data on the spectral hard state of black-hole binaries show the 
high-energy cut-off around ~ 100 keV (e.g., Grove et al. 1998; 
Wardzinski et al. 2002; Zdziarski & Gieiiinski 2004; Done, 
Gierlinski & Kubota 2007; Joinet, Kalemci & Senziani 2008). 
When this cut-off is parameterized by the maximum in the EdF/dE 
plots, it is in the range of ~ 50-200 keV (for spectra with the photon 
index < 2, which is the case in the hard state). 

We compare these results with our model spectrum, see Fig. 
[3{d). The maximum in the EdL/dE spectrum of the initial, not self- 
consistent, solution is at ^ 700 keV if GR effects are neglected 
and Comptonization is local (the long-dashed curve), which then 
moves to 600 keV when GR effects are taken into account (the 
short-dashed curve). The peak energy is at ^ 500 keV for the self- 
consistent solution when both global Comptonization and GR ef- 
fects are neglected (the dot-dashed curve). The value of this peak 
energy is then very similar to the value of the electron temperature, 
kT s ~ 700 keV (Fig.|3^) at the radius of the peak dissipation of the 
corresponding model (Fig. [3]:). The final dynamic solution has the 



dissipation rate much more centrally concentrated than the initial 
one, see Fig.[3jc), and then the GR effects reduce the peak energy 
much more, to a 300 keV (the solid curve). Thus, taking into ac- 
count Comptonization self-consistently, i.e., as the global process 
in the entire hot disc, and including the GR effects reduces the peak 
energy of the spectrum by a factor of ^ 2/3. This is because the 
global Comptonization increases the spectral contribution of outer 
flow regions while the GR effects reduce the contribution of in- 
ner regions. The GR effects here include three parts (see discussion 
in Section [5J, namely the gravitational and Doppler redshift, and 
the preferential capture of most energetic photons (enhanced by 
the kinematic collimation) by the black hole. The last effect is the 
dominant one. 

However, this decrease is still not sufficient to bring the theo- 
retical cut-off down to the observed values. This problem appears 
to be shared by the hot accretion flow model in general, see, e.g., 
Yuan & Zdziarski (2004). As stated above, the cut-off energy before 
correcting for the GR effects approximately equals the temperature 
of the flow in the region where most of the radiation is emitted. 
Emission of this region can then be roughly approximated as that 
of a one-zone thermal-Comptonization model. Fits of this model 
to hard-state spectra of black-hole binaries and of Seyfert galax- 
ies show kT s ^ 50-100 keV, see, e.g., the compilation in Yuan & 
Zdziarski (2004). Approximately, kT c of the one-zone model equals 
the energy of the maximum of its EdE/dE spectrum. The other pa- 
rameter of this model is the Thomson optical depth, r-r. For a given 
X-ray spectral index, the Compton parameter, y = 4T e T T , is ap- 
proximately constant. Also, kT e is anticorrelated with tt for a given 
local dissipation rate in a hot flow, reflecting the varying power per 
electron. Fits of the one-zone thermal-Compton model yield t t ~ 1 
(e.g., Yuan & Zdziarski 2004) whereas the vertical optical depths 
of the flow are tj <k 1, see Fig.|3tb), which expresses the same dis- 
crepancy as r c of the flow being too high. As stated above, taking 
into account the GR effects in the global-Compton model improves 
the agreement significantly (pointing to the inadequacy of Comp- 
tonization models neglecting GR), but not sufficiently. 

On the face of it, this discrepancy might indicate that the hot 
flow model is not applicable to the hard state of black hole binaries 
and Seyfert galaxies. However, we note that the present calcula- 
tions assume a very strong outflow, equation (|3) with s = 0.3, mo- 
tivated by the modelling of the Galactic Centre (Yuan et al. 2003). 
On the other hand, our assumed dimensionless accretion rate is » 
that of the Galactic Centre source, whereas the relative strength 
of outflows appears to decrease with the increasing accretion rate 
(Gallo, Fender & Pooley 2003; Fender, Belloni & Gallo 2004), as 
also indicated by a relatively small change of the X-ray bolomet- 
ric luminosity during hard/soft state transitions (e.g., Zdziarski et 
al. 2004). This is plausible from a theoretical point of view because 
the Bernoulli parameter becomes smaller with the increasing accre- 
tion rates (Yuan, Cui & Narayan 2005; Bu & Yuan in preparation). 
Thus, it is likely that the fractional strength of the outflow, 82 per 
cent for our assumed s, R ollt and M (with 18 per cent of M cross- 
ing the horizon), may in reality be much weaker. Setting 5 = 
leads to t t being a few times higher and T e about a factor of two 
lower, see Fig. [2] Although the calculations in Fig.|2]are done with 
the local analytical treatment of Comptonization, and are thus not 
self-consistent, they correctly predict the direction of the changes 
of the flow parameter. If the characteristic T e of the self-consistent 
flow were reduced by the same factor of two with respect the model 
with s = 0.3, the peak of the EdL/dE spectrum would also go down 
by a similar factor, bringing it to the observed range and resolving 
the discrepancy with the data. 
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Our second assumption has been of a strong viscous heating 
of electrons, S = 0.5. As seen in Fig. [2] this has a relatively minor 
effect on T c and t t . We also note that the half-depth Thomson opti- 
cal depth of a slab is the quantity closest to that of the half-depth of 
the flow, and should preferably be used when comparing accretion 
flow models with one-zone thermal Comptonization models. For a 
given X-ray spectral index, the optical depth is somewhat lower for 
a slab geometry than for spherical one. 

We note that there is a number of additional effects that can 
further reduce T e and increase t t , possibly allowing the hot flow 
model to be in agreement with the full range of the observed hard- 
state spectra. First, an increase of the black-hole spin reduces the 
radial velocity and increases the density of the flow. The effect is 
rather strong, as can be seen, e.g., in fig. 5 of Gammie & Popham 
(1998). Global Comptonization in the Kerr metric will be studied 
in detail in our forthcoming work. Second, the presence of mod- 
erate large-scale toroidal magnetic fields in the accretion flow sig- 
nificantly reduces T e , as shown by Bu, Yuan & Xie (2009). Third, 
electron cooling will also be significantly enhanced in either a two- 
phase flow, with cold clouds mixed with the hot flow, or, fourth, 
at presence of an inner collapsed disc, with both effects happening 
above some critical accretion rate (Section[5] see also Y01, Yuan 
2003). 

Furthermore, the cut-off energy in the hard state of black- 
hole binaries is observed to decrease with the increasing luminos- 
ity (Wardzihski et al. 2002; Yamaoka et al. 2006; Yuan et al. 2007; 
Miyakawa et al. 20080 Also, the hot-flow characteristic temper- 
ature goes down with the increasing accretion rate, in agreement 
with the observations. In our case, the bolometric luminosity is 
only 0.007Le (Section[5j, and thus our model spectrum should be 
compared with the hard state at the correspondingly low L (which 
is somewhat below, e.g., the hard-state range of L of Cyg X-l of 
- 0.01-0.02Le, Zdziarski et al. 2002), which spectra are likely 
to have the cut-off energies > 100 keV. This would be in agree- 
ment with our results, after the correction for the outflow discussed 
above. 

Malzac & Belmont (2009) have also pointed out that the ion 
temperatures implied by observations of the hard state are much 
lower than those typical for ADAF models. However, the ratio of 
TJT t , calculated by Malzac & Belmont (2009) directly from the 
formula for Coulomb energy transfer from ions to electrons, is ap- 
proximately oc Tj 2 , and the discrepancy pointed out by those au- 
thors occurs at tt ~ 1. If we take instead the low values of tt 
obtained in ADAF models, there is an agreement between their es- 
timate and the hot flow models. 

On the other hand, we stress that using equation $3$ with j- = 
0.3 does not weaken our conclusions related to the role of global 
Comptonization. Instead, if the outflow is weaker for a given M , 
the radiation generated at small radii and subsequently received at 
large radii will be stronger, thus the global Comptonization effect 
will become even more important than for our chosen assumption. 
This would further reduce the average photon energy of the emitted 
spectrum. 

Another related issue is our determination of the value of the 
critical accretion rate, at which an inner part of the hot flow col- 
lapses. We find it corresponds to a relatively low L ^ O.OILe. Lu- 



Note that temperatures as low as ~ 20 keV were claimed in the one-zone 
model fits in Miyakawa et al. (2008). However, they were obtained using 
the non-relativistic Comptonization model of Sunyaev & Titarchuk 1980), 
which is not valid for spectra extending above ~ 50 keV. 



minosities of the hard state of black-hole binaries are commonly 
above this L (e.g, Done et al. 2007; Zdziarski et al. 2002). Again, 
our determination is for the assumed strong outflow, and relaxing 
this assumption may increase that critical luminosity. On the other 
hand, the presence of a collapsed geometrically thin disc close to 
the horizon would explain the finding of relativistically broaden Fe 
Kor line in the hard state of black-hole binaries (e.g., Miller et al. 
2006; Miller 2007; but see Done & Diaz Trigo 2009 for a critical 
view), whose origin would have otherwise been in conflict with the 
hot accretion flow model of the hard state. 

Finally, we note that alternative models for the hard state have 
been proposed. A very interesting recent model is non-thermal, in 
which the power supplied to electrons (with the optical depth of 
Tt ^ 1) goes into their acceleration into a power-law distribution 
(Poutanen & Vurm 2009; Malzac & Belmont 2009). Then, syn- 
chrotron self-absorption and Coulomb interactions efficiently ther- 
malize the electrons provided any blackbody emission is weak. 
This radiative one-zone model fits the hard-state data, e.g., of Cyg 
X-l, very well. An important issue here is the location of the non- 
thermal plasma. It cannot be a corona as the model constrains any 
disc blackbody emission irradiating the plasma to be very weak. If 
it is an inflow, it would require an efficient mechanism of electron 
acceleration (instead of heating) with a large S, and a low radial ve- 
locity to achieve the required large t t . It appears unclear whether 
such a model is possible dynamically. 



7 CONCLUSIONS 

Hot accretion flows, both ADAF and LHAF, are optically thin in 
the radial direction, and thus the global Compton scattering is im- 
portant. Its effect on the dynamics of the accretion flow has been in- 
vestigated previously by an analytical approach in YXO09. In this 
paper we revisit this problem using the more accurate MC simula- 
tions, and focusing on the inner region of the flow. We confirm that 
the global Compton scattering effect is dynamically very impor- 
tant. We obtain the final self-consistent solution using the iteration 
approach. Our main results, in particular showing the differences 
between the initial and final iterations, are presented in Fig. [3] The 
global GR cooling rate of the final self-consistent solution is more 
centrally concentrated than that of the initial solution, but less cen- 
trally concentrated than the local cooling rate of the final solution. 

To evaluate the effect of global Compton scattering, we also 
compare the global GR Compton scattering model with two local 
models, in the sphere and slab geometries. In these models, only the 
local scattering is taken into account. For both the initial and final 
self-consistent solutions, we find that none of the local models can 
adequately reproduce the global model, although the slab model is 
significantly better than the spherical one. This is due to the fact 
that the radial optical depth of the hot flow is much higher than the 
vertical one. 

Specifically, the local slab model underestimates the cooling 
rate in an outer region, > 200R S , by up to two orders of magnitude 
with respect to the global GR MC results, while it overestimates the 
cooling rate by up to factor of a few in an inner region, < 20R g . This 
differences occur because of the transfer of seed photons between 
different radii, mostly from small to large radii. We also find that 
the main dynamical effect in our self-consistent treatment is the 
global Comptonization, whereas the GR and bulk motion effects 
are dynamically less important. 

We also compare the spectra of the initial and final solutions 
in both the local and global models. We find that the shape of the 
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spectrum of the final self-consistent solution changes only slightly 
(mostly in the reduced high-energy cut-off) compared to that of the 
initial solution, but the normalization decreases by a factor of ~ 2. 
We also find that the spectral differences between global and local 
slab models are much less significant than those in the radial profile 
of the cooling rate. This is due to the averaging of the local spectra 
over all radii. 

We have compared our results to the spectra of the hard state of 
black-hole binaries, see Section[6] We find an overall rough agree- 
ment, though our models have the characteristic electron temper- 
atures and the characteristic Thomson optical depths higher and 
lower, respectively, than those required by the observational data. 
However, this discrepancy can, most likely, be resolved if the disc 
outflow in the luminous hard states is much weaker than that as- 
sumed by us. Additional factors that would allow the model be in 
agreement with the full range of the observed spectra are the spin of 
the black hole, large-scale magnetic fields, and the presence of an 
additional electron cooling mechanism, e.g., soft photons from cold 
clouds embedded in the hot flow, or from an inner collapsed thin 
disc. We stress that global Comptonization is crucial in bringing 
the hot accretion flow model to agreement with the data. This is due 
to the self-consistent model having the dissipation much more cen- 
trally concentrated than that of local models, which in turn causes 
the GR effects to reduce the spectral cut-off energy much more for 
the former than the latter. 
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